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Abstract 

We solve the generalized relativistic harmonic oscillator in 1+1 dimensions, i.e., including a linear 
pseudoscalar potential and quadratic scalar and vector potentials which have equal or opposite 
signs. We consider positive and negative quadratic potentials and discuss in detail their bound- 
state solutions for fermions and antifermions. The main features of these bound states are the same 
as the ones of the generalized three-dimensional relativistic harmonic oscillator bound states. The 
solutions found for zero pseudoscalar potential are related to the spin and pseudospin symmetry 
of the Dirac equation in 3+1 dimensions. We show how the charge conjugation and 7^ chiral 
transformations relate the several spectra obtained and find that for massless particles the spin 
and pseudospin symmetry related problems have the same spectrum, but different spinor solutions. 
Finally, we establish a relation of the solutions found with single-particle states of nuclei described 
by relativistic mean-field theories with scalar, vector and isoscalar tensor interactions and discuss 
the conditions in which one may have both nucleon and antinucleon bound states. 

PACS numbers: 21.10.Hw, 21.60.Cs, 03.65.Pm 



I. INTRODUCTION 

Recently, there has been a wide interest in relativistic potentials involving mixtures of 
vector and scalar potentials with opposite signs. The interest lies on attempts to explain 
the pseudospin symmetry in nuclear physics. Chen et al. [l|, using a Dirac Hamiltonian 
with scalar Vs and vector Vy potentials quadratic in space coordinates, found a harmonic- 
oscillator-like second order equation which can be solved analytically for A = Vf — Vs = 0, as 

n 

considered before by Kukulin et al. '2\, and also for S = Vi+K = 0. Very recently, Ginocchio 
solved the triaxial, axial, and spherical harmonic oscillators for the case A = and applied it 
to the study of antinucleons embedded in nuclei [3[ . The case S = is particularly relevant 
in nuclear physics, since it is usually pointed out as a necessary condition for occurrence of 
pseudospin symmetry in nuclei |^-|8|. Also the observed isospin dependence of pseudospin 
can be explained by the effect of the p potential on the shape of the S potential |9|-[l3|. As 
shown by Bell and Ruegg 1J|, the A = and S = cases correspond to a SU{2) symmetry 
of the Dirac Hamiltonian with only scalar and vector potentials which is independent of the 
particular shapes of these potentials. 

A generalized harmonic oscillator (HO) for spin 1/2 particles that includes not only those 
combinations of scalar and vector potentials, but also a linear tensor potential, giving rise 
to the so-called Dirac oscillator, has also been considered in ,15]. The research of this kind 
of interaction was started by Ito 'ig] and has been revived by Moshinsky and Szczepaniak 
|l2| (see |l5| for a complete reference list). 

In Ref. [l5| a special attention was paid to the case when the scalar and vector potentials 
are such that S = 0, which, as stated before, has been related to the nuclear pseudospin 
symmetry. There it has been concluded that only when the tensor potential is turned 
on one gets negative energy bound-state solutions along with the usual positive energy 
solutions. Those negative energy solutions are important to describe the possible existence 
of antinucleons in nuclei JS|. When the tensor potential is absent the special conditions 
relating the scalar and vector potentials (Vs = ±Vt) are needed to have a standard positive 
HO potential which is only able to bind fermions, and therefore exclude the negative bound- 
state solutions from the spectra. Actually, as was discussed in another context in [l8[, 
with an appropriate choice of signs, potentials fulfilling the relations Vg = ±\4 are able to 
bind either fermions or antifermions. As an example, we will see later that a (negative) 



harmonic potential, with no tensor potential, binds antiferniions and does not allow positive 
bound-state solutions. 

Closely related to this is the fact that, in the nucleus, the charge-conjugation transfor- 
mation relates the spin symmetry of the negative bound-state solutions (antinucleons) to 
the pseudospin symmetry of the positive bound-state solutions (nucleons) ,19]. This has 
also been discussed recently in j3j], analyzing the HO for antinucleons with spin symmetry 
(A = 0) and without the tensor interaction. Again, the charge conjugation connects this 
problem with the HO for nucleons with pseudospin symmetry (S = 0), but in this case the 
positive harmonic binding potential is A. Note that, because the harmonic potential does 
not vanish for large distances but is a confining potential, it is possible to obtain particle 
bound-state solutions in the exact pseudospin limit 15]. Therefore, we believe that this 
connection between spin and pseudospin symmetry obtained by charge conjugation deserves 
to be more explored. In addition, after the pseudospin symmetry has been revealed as a 
relativistic symmetry by Ginocchio |J|, a link of it to chiral symmetry has been suggested 
based on sum rules calculations in relativistic chiral models that support a small value for 
E J20]. However, until now, a clear connection of pseudospin and chiral symmetries has not 
been established. 

The main motivation of this paper is to shed some light in the relation between spin 
and pseudospin symmetries by means of charge-conjugation and 7^ chiral transformations. 
In order to do so, we solve the mixed scalar-vector-pseudoscalar HO potential in 1-1-1 di- 
mensions, taking advantage of the simplicity of the lowest dimensionality of the space-time. 
This approach is equivalent to consider fermions in 3-1-1 dimensions which are restricted to 



move in one direction J2l|. We explore the spectra when it is possible to obtain analytical 
solutions, i.e., in the particular cases when A = and S = 0. As referred before, if the 
pseudocalar term is turned off, these cases correspond, respectively, to have exact spin and 
pseudospin symmetry in 3-1-1 dimensions. We explore all the possible signs of the quadratic 
(A or S) and linear (pseudoscalar) potentials, thus paying attention to bound states of 
fermions and antiferniions as well. We compare both cases A = and S = to establish the 
charge-conjugation connection discussed above in the presence of the pseudoscalar term. We 
also consider the case of zero mass to look for the connection between spin and pseudospin 
symmetries by means of the chiral transformation. 

An added motivation for this paper is given by the recent demonstration that the rela- 



tivistic 3+1 HO with scalar and vector potentials can describe successfully the heavy nucleus 



spectrum 2^. The parameters of the HO are determined by fitting the scalar and vector 
potentials derived from Relativistic Mean Field nuclear calculations (RMF). The ^°^Pb neu- 
tron single-particle levels obtained in the HO are very similar to the RMF results. The 
pseudospin symmetry is shown to be almost satisfied for a heavy nucleus in the HO descrip- 
tion. The possibility of breaking perturbatively this symmetry in the 3+1 HO by a tensor 



interaction (in our 1+1 case a pseudoscalar coupling) has been discussed in J23| and can be 



included in a calculation like the one of Ref. 



22| in order to improve the results. Thus, the 



study of all the possible eigenenergies of the 1+1 HO, presented in this paper, considering 
not only the positive energy solutions already obtained for the 3+1 case in 15| but also the 
negative ones, can be applied to analyze the existence of antinucleons in nuclei and their 
spectrum. 

The particle and antiparticle spectra we obtain for the cases E = and A = demon- 
strate the charge-conjugation connection referred above. They also show that the spectra 
of massless A = states are degenerate with massless S = states. Further imposing chi- 
ral symmetry to these states, such that they would become eigenstates of 7^, would imply 
that all potentials would be zero and thus one would have just massless free fermions. Our 
conclusions are valid in the case of 3+1 dimensions, because, as will be shown later, one 
basically has to add the angular momentum quantum numbers £ and n to the coefficients of 
the functions appearing in the 1+1 eigenvalue equations in order to get the corresponding 
3+1 equations. 

This paper is organized as follows. In Sec. |H] we present the general Dirac equation in 
1+1 dimensions with a potential with a completely general Lorentz structure. The effect of 
charge-conjugation and the 7^ chiral transformations upon the general 1+1 Dirac Hamilto- 
nian is discussed in Sees. IIIAI and IIIBI respectively. In Sees. Ill CI and IIIDI we present the 
equations of motion and discuss the isolated solutions, i.e., the solutions out of the Sturm- 
Liouville problem, for A = and S = 0. The nonrelativistic limits of the A = and S = 



cases are discussed in Sec. IIIE| where we show that when S = there is no nonrelativistic 
limit when the pseudoscalar potential vanishes or is very small. In Sec. Illll we introduce the 
general relativistic harmonic oscillator in 1+1 dimensions for A = by assigning a quadratic 
potential kix'^/2 for S and a linear potential k2X for the pseudoscalar potential Vp. There 
we also present the eigenvalue equation and the wave functions for bound states, discuss 



the Dirac oscillator case (A = E = and Vp ^ 0) and show how one can use the charge- 
conjugation and chiral transformations to get the eigenenergies for the S = case from the 
A = case. The following section is devoted to a detailed analysis of the solutions of the 
A = relativistic HO for all signs of the parameters ki and k2. Finally in Sec. |3we discuss 
the spectrum for S = and massless fermions and in Sec. I VII we draw the conclusions. 

II. THE DIRAC EQUATION IN 1+1 DIMENSIONS 

The two-dimensional Dirac equation can be obtained from the four-dimensional one with 
the mixture of spherically symmetric scalar, vector and anomalous magnetic-like (tensor) 
interactions. If we limit the fermion motion to the x-direction {py = Pz = 0) the four- 
dimensional Dirac equation decomposes into two equivalent two-dimensional equations with 
2-component spinors and 2x2 matrices i21||. Then, there results that the scalar and vector in- 
teractions preserve their Lorentz structures whereas the anomalous magnetic-like interaction 
turns out to be a pseudoscalar interaction. Furthermore, in the 1+1 world there is no an- 
gular momentum so that the spin is absent. Therefore, the 1+1 dimensional Dirac equation 
allows us to explore the physical consequences of the negative-energy states in a mathe- 
matically simpler and more physically transparent way. In this spirit the two-dimensional 
version of the anomalous magnetic-like interaction linear in the space coordinate has also 
received attention [2^- [2^. Later this system was shown to be a particular case of a more 
general class of exactly solvable problems f30| . 

Let us begin by presenting the Dirac equation in 1+1 dimensions. The Dirac equation 
for a particle with mass m and a general potential V, to be defined later, is formally the 
same as its 3+1 counterpart: 

ihcYd^,^ - mc^^ - 7°V^ = . (1) 

Here /i = 0, 1 and the 2x2 matrices 7^^ obey the usual relations {7'^,7'^} = 2g'^'^I, where 
/ is the 2x2 identity matrix. The metric tensor and the covariant derivative operator are, 
respectively, g^'^ = diag(l,— 1) and 9^ = {d/d{ct) , d/dx). The Hamiltonian is defined in 
the usual way, such that Eq. (P) is written as 

ih^— = H-^^ = {cap + dmc^ + V)'^ , (2) 

ot 



where p = —ihd/dx. The positive definite function |\l/p = \1/''^\1/, satisfying a continuity 
equation, is interpreted as a position probabihty density and its norm is a constant of 
motion. This interpretation is completely satisfactory for single-particle states J31|. The 
traceless matrices a and (3 are defined by a = 7^7^ and /3 = 7°, and obey the relations 
a'^ = f3^ = /, {a, 13} = 0. We set V to be 

V = IVt + l3Vs + aVsp - 1(3^% , (3) 

where 7^ = 7^*7^ = a. This is the most general combination of Lorentz structures because 
there are only four linearly independent 2x2 matrices. The subscripts for the terms of the 
potential denote their properties under a Lorentz transformation: t and sp for the time and 
space components of the 2- vector potential, s and p for the scalar and pseudoscalar terms, 
respectively. 

If the terms in the potential V are time- independent, the Dirac equation becomes 

Hi^ = E^ (4) 

where 

H = cap + pmc^ + I Vt{x) + aVsp{x) + pV^ix) - ip-f%{x) , (5) 

^(x) = e^^*^(a;,t) (6) 

and E is the fermion energy. To have an explicit expression for the a and (3 matrices one 
can choose 2x2 Pauli matrices which satisfy the same algebra. We use (3 = a^, a = ai and 
thus P'j^ = ia2- 

The Hamiltonian (0) is invariant under the parity operation, i.e., when x -^ —x, Vsp{x) 
and Vp{x) change sign, whereas Vt{x) and Vs{x) remain the same. This is because the 
parity operator is P = exp(iry)PoC"3, where r^ is a constant phase and Pq changes x into 
—X. Since this unitary operator anticommutes with a and jSj^, they change sign under a 
parity transformation, whereas / and (3, which commute with P, remain the same. When 
one writes down the explicit equations of motion in terms of the components of the spinor 
ip, the combinations Tj = Vt + Vs and A = Vt — Vs oi the vector and scalar components 
arise naturally. Therefore, it is convenient to rewrite the Hamiltonian (0) in terms of these 
potentials. We have 

H = cap + pmc" + aV^p + ^^S + ^^-^ A - 1/3^% . (7) 



A. Charge conjugation 

The charge-conjugation operation changes the sign of the electromagnetic interaction, 
i.e., changes the sign of the time and space potentials in (j2)). This is accomplished by the 
transformation (see, for example, Itzykson and Zuber |32l |) 

^ — > ^, = C^^ = C7°'^^* , (8) 

where T denotes matrix transposition and C is a matrix such that C^y'^'^C"^ = —7'^. In 
1 + 1 dimensions one matrix that satisfies this relation is 

C = e*^a/? . 

If we choose the phase factor e*^ equal to 1, we have 

^^ = C^T = a^* . (9) 

After applying this charge-conjugation operation to the Dirac equation Q, the time- 
independent Dirac equation becomes 

H,i,, = -Eij, , (10) 

where ipc = oitp* and He is given by 

H, = cap + (3m^ -IVt- aVsp + pVs + i(3-i% . (11) 

In terms of the potentials A and S, this Hamiltonian reads 

H, = cap + /3mc^ - aVsp - ^^-^A - ^^-^S + if3y% . (12) 

One sees that the charge-conjugation operation changes the sign of the energy and of the 
potentials Vt, Vsp and Vp. In turn, this means that S turns into —A and A into — S. 
Therefore, to be invariant under charge conjugation, the Hamiltonian must contain only a 
scalar potential. 

B. Chiral transformation 

The chiral operator for a Dirac spinor is the matrix 7^, and we will call "chiral trans- 
formation" the transformation associated with it. Thus, the transformed spinor is given by 

8 



\1/^ = 7^\1/ and the transformed Hamiltonian H^ = 'j^H'-f^. Since 7^ anticommutes with /3, 
the time-independent chiral-transformed Dirac equation is 

H^^x = Ei^x , (13) 

where H^ is given by 

H^ = cap- pmc^ + IVt + aV^p - /?K + ipj% (14) 

or 

H^ = cap- (3mc^ + aVsp + ^^-^A + ^-^^ + ^/^T^p (15) 

in terms of S and A. This means that the chiral transformation changes the sign of the 
mass and of the scalar and pseudoscalar potentials, thus turning S into A and vice-versa. 
A chiral-invariant Hamiltonian needs to have zero mass and Vg and Vp zero everywhere. 

C. Equations of motion 

The space component of the 2-vector potential in Q can be absorbed into the wave 
function by defining a new spinor if) such that 

^ = e-^^V , (16) 

in which A = (l/hc) J^ Vsp{x') dx', since we have Hip = e~*^(// — aVsp)ip. From this point 
on, we will refer to Vt as simply a vector potential, following the common usage of this term 
(usually denoted by Vy). If we now write the spinor ip in terms of its components. 



^ = , (17) 




the Dirac equation gives rise to two coupled first-order equations for ip+ and ■?/'_: 

- ihcipi + mc^ilj+ + S^+ - i Vpip- = EiIj+ (18) 

-ihcifj'j^ - mc^%l)_ + /\%l)_ + i Vpilj+ = Eip_ , (19) 

where the prime denotes differentiation with respect to x. In terms of tp^ and ^_ the 
spinor is normalized as /_ dx (IV'+P + IV'-P) = 1) so that ip^ and ip- are square integrable 



functions. It is clear from the pair of coupled first-order differential equations ()18p and ()19|) 
that ?/'+ and tp- have opposite parities if the Dirac equation is covariant under x —>■ —x. 
Under the charge-conjugation and chiral transformations the spinor (|17|) becomes 



(20) 




and 




(21) 



respectively. 



D. The Sturm-Liouville problem and isolated solutions 

Using the expression for ijj_ obtained from (fTU|) with E ^ —mc^ + A, viz. 

i,_ = _,M.ZJ% (22) 

and inserting it in ()18j) one arrives at the following second-order differential equation for ?/;+: 
-h\^i;l+hc^''^^^i^±^P^+ \V^ + hcV' -{E- mc^ - S)(E + mc^ - A)l ^+ = . (23) 
In a similar way, using the expression for ip^ obtained from (fTH|) with E ^ mc^ + S, viz. 

and inserting it in pfjj) one arrives at the following second-order differential equation for ip^: 

V, lb- + hcib' 

- h^^ip'L - hcE'-^ ^— ^ + \V^ - hcV' -(E- mc^ - S)(E + mc^ - A)l ^_ = (25) 

E — mc' — E '- ^ ^ 

For A = with E ^ -mc"^, ^ and (j231) reduce to 

*-=-.^^le^, (26) 

E + mc^ 



2 



2m + 



E + mc^ ^ Kf hV' 



S + -^ + 



p I p 



2mc^ 21710^ 2mc 

and for S = with E ^ mc^, 1^^ and ^ reduce to 
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2m 



A+ ^ 



hv; 



2^4 



^_ 



E^ - m^c 
2mc^ 



■^. 



(29) 



2mc^ 2mc^ 2mc 

Either for A = with E ^ —mc^ or S = with E ^ m(? the solution of the relativistic 
problem is mapped into a Sturm-Liouville problem in such a way that solution can be found 
by solving a Schrodinger-like problem. If one considers potentials Vt and K quadratic in x 
with Vt = ±Vs and a potential Vp linear in x one obtains Schr'odinger-like equations for the 
harmonic oscillator potential. 

The solutions for A = with E = —mc^ and S = with E = mc^, excluded from the 
Sturm-Liouville problem, can be obtained directly from the original first-order equations 
(UHl) and (Uni). They are 



^p+ = ijj\ exp 



he 



+ 



dy- 



he 



> , for A = with E 



he 



(S + 2me^) V'^ 



-me 



(30) 



and 



^_ 



/(O) 

tp_ exp 



- / dy 



do) 



,(0) 



Vp{y) 

he 



-— (A-2mc^)^_ 
he ^ ' 



> , for S = with E 



mte 



(31) 



where ■?/'+ and ?/'_ are normalization constants. Since, by parity conservation, V^ must be 
an odd function of x, its integral must be an even function of x. Thus, if the integral of V^ 
has a definite sign, its overall sign can determine whether the isolated solutions ()3U|) and 
fl31|) are bound states. In particular, if Vy and S or A are powers of x, one sees from each 
pair of equations (jnUj) and (j^ that, in general, one of the spinor components has to be zero, 
either ?/'+ in the first case (sign of the V^ coefficient positive) or ?/'_ in the last case (sign of 
the Vp coefficient negative). This is the case of the harmonic oscillator potentials, as we will 
see in the next section. However, it is possible find solutions with the Vp coefficient negative 
in the first case (A = 0, -E = — ?7ic^) and with the Vp coefficient positive in the second case 
(S = 0, -E = m(?\ provided a certain relation between the coefficients of all the potentials 
is satisfied, as we will see later. 
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E. The nonrelativistic limit 



In the nonrelativistic limit, E ~ mc^, Eq. ()23|) becomes 



<-(2^;^,'''"'--h-^^'i'"-^&-"^ 



+ SV^+ = £^+ , (32) 



2mc2 - A ^ 

where S = E — mc^. Since, in the nonrelativistic regime, for the range of values of x in 
which the wave function is not negligible, |A(x)| <^ mc^, from Eq. (f^^ we see that ip- is of 
order v/c <^ 1 relative to ip+, provided, of course, one has also |V^(a;)| <^ mc^. On the other 
hand, since in these conditions 2r?7,c^ — A ~ 2mc^ and the A' term in Eq. (jH^ is suppressed 
relative to the other terms, we also see that 4'+ obeys a Schrodinger equation with binding 
energy equal to S = E — mc^. 

We shall consider now the nonrelativistic limit in the special cases A = and S = 0, 
which we will also consider later in the paper when we solve the Dirac equation for harmonic 
oscillator potentials. 

When A = 0, Eq. §^ becomes 

Here we clearly see that S plays the role of a binding potential in the nonrelativistic limit 
and that Vp gives rise to effective binding potentials proportional to V' and Vp. This means 
that even a pseudoscalar potential unbounded from below could be a confining potential. 
Note that if Vp/lmc^) is of the same order that Ti/lmc^) in a l/(mc^) expansion within 
the classically accessible region (positive kinetic energy) where the wave function is not 
negligible, and if its variation over this region is such that h/ {2'mc)Vp ~ {Vp changes very 
little over a distance of the order of a Compton wavelength), Vp gets suppressed relative to 
E. 

For the case E = and if |A|/(mc^) is very small in the classically accessible region, 
Eq. (jn21) becomes 

If V^ = or if, as explained before, the Vp and Vp terms are higher-order terms in a 
l/(mc^) expansion within the classically accessible region, Eq. (jH^ becomes a free particle 
Schrodinger equation. This is in agreement with what was found in Ref. |l5|, namely that 
when T, = Vp = and A is a three-dimensional harmonic oscillator potential, there is no 
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nonrelativistic limit. In this case we were able to show that this is true in a more general 
framework. A model for a relativistic particle in which the vector and scalar potentials 
are such that S = 0, and where the pseudoscalar potential is suppressed, is intrinsically 
relativistic, as far as bound states are concerned. 

This last result deserves some more comments. The fact that in a Dirac equation one 
can have Lorentz scalar and pseudoscalar potentials leads to remarkable results that are 
at odds with what is known from nonrelativistic quantum mechanics, where the Lorentz 
structure plays no role, since one has only potentials which couple to the energy, i.e., behave 
as time components of a relativistic vector. Among these results is the fact that the Dirac 
equation is not invariant under a simultaneous shift of the energy, the scalar potential. 



and the pseudoscalar potential. It has &. 
)seudoscalar screened Coulomb potential 



ready been verified that a constant added to the 

-I 

33| or to the pseudoscalar inversely linear potential 

3J| is physically relevant and plays a crucial role in ensuring the existence of bound states. 
As an example of another notable result, it is well known that a confining potential for 
a nonrelativistic fermion cannot confine it in the relativistic regime when the potential is 
considered as a Lorentz vector (see, e.g.., pill), simply because there is pair creation and the 
single-particle picture does not hold anymore. It is surprising that the converse is also true, 
i.e., that relativistic binding potentials may not bind in the nonrelativistic limit. This is 
the case oiVt = —Vg (S = 0) referred before, and is basically due to the fact that vector 
potentials couple in a different way than scalar potentials in the Dirac equation, whereas 
there is no such distinction in the Schrodinger equation. Actually, as we have seen before, E, 
which is the sum of a vector and a scalar potential, plays the role of a binding potential in the 
nonrelativistic limit, whereas a weak enough A cannot bind a fermion with nonrelativistic 
energy if S = 0. Note that for potentials going to zero at infinity, like mean-field nuclear 
potentials, there are no bound states when E = 0, even if A is not small compared with 
fermion rest energy (4|-(lll|. 

Regarding antiparticles, an attractive vector potential for a particle is, of course, repulsive 
for its corresponding antiparticle, and vice versa. However, the attractive (repulsive) scalar 
potential for particles is also attractive (repulsive) for antiparticles. This is expressed in 
the change A -^ — S and S -^ —A under charge conjugation, as referred in Sec. IIIAI This 
means that an antiparticle with nonrelativistic energies {E ^ —mc^) will not be bound when 
A = 0, provided Vp is small enough. 
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III. THE RELATIVISTIC HARMONIC OSCILLATOR 

Let us consider 

E = -kix^, A = 0, Vp = k2X (35) 

As we have seen in the last section, the chiral transformation performs the changes 
A — > S, S — i> A, ?7i — ^ —m and Vp — > —Vp. Morever, since 7^ interchanges the upper 
and lower components (see ((23), the resulting pair of transformed equations of motion are 
formally the same, so that their solutions have the same energy eigenvalues. This symmetry 
can be clearly seen from the two equation pairs (J22j) - (j23j) and (I24|) - (j25|) as well as from 
the isolated solutions (jHn|l and (pTT| which are converted into each other by this kind of 
transformation. This means that one can restrict the discussion to the E case (A = 0). The 
results for the case when A = ^kix^, S = 0, Vp = k2X can be obtained immediately by just 
changing the sign of m and of k2 in the relevant expressions. 

The Dirac spinor corresponding to the isolated solution with E = —mc^ is obtained 
from (jHOI)- Only for k2 > there is a normalizable Dirac spinor and, as commented 
before, the upper component vanishes whereas the lower component is given by ip- = 
■ip_ exp{—k2x'^/{2hc)), regardless the value of ^2. For E ^ —mc^, Eq. P7j) takes the 
form 

- ^^: + Ikx'^P+ = ^^+ (36) 



where 



' E + mc^ 



K=^{—^ h + H] (37) 



~ E"^ — rn^c^ hk2 .„„. 

E = . 38) 

2mc^ 2mc 

The well-behaved solution for ()36|). with K necessarily real and positive, is the well-known 

solution of the Schrodinger equation for the nonrelativistic harmonic oscillator (see, e.g., 
Q): 

2P+ = N^H^{Xx)e-^"^"/' (39) 



E-h + l)N^ (40) 



14 



where n = 0,1,2,..., Nn is a normalization constant, Hn (Ax) is a n-th degree Hermite 
polynomial and 

The lower component of the Dirac spinor is obtained from ()26|) . When one uses the 
recursion relations for the Hermite polynomials (see, e.g., J2S|) one gets 



ihcXNr, 



H^+^ [Xx) - ( 1 - -^ ) (\x) H^ {\x) 



e-^^ /2 (42) 



E + mc^ 

Since the Hermite polynomial if„ (Ax) has n distinct zeros one may conclude that ■?/'+ has 
n nodes, and the expression for ip^ suggests that it has nil nodes, depending on the sign 
^2. In fact, if fci = 0, this last expression reduces to 



ihcXNr, 

ip. 



i7„+i(Ax)-(l-|^)(Ax)i7„(Ax) 



e-"'^-'/2 , (43) 



E + mc^ 

showing that ?/'_ is proportional to Hn+i{\x) if /c2 is positive. If ^2 < 0, the recursion 
relations of Hermite polynomials imply that i})-. is proportional to Hn~i{\x) (see J36|). 

From Eqs. (|H7j). (I^Hj) and ^^ we can get the quantization condition for the Dirac 
eigenenergies: 

E^ - m^c^ = {2n + 1) he W :^±^A;i + kl + hc ^2 (44) 

The nonrelativistic limit is reached when h^\ki\/{nr'c^) ^ 1 and h\k2\/{m?c') ^ 1. In this 
case E ~ mc^ for small quantum numbers, and Eq. (j44|) becomes 



so that ki is restricted to ki > —k'^/lmc^). 

In general there is no requirements on the signs of ki and k2, except that, as stated before, 
i^ > 0, and therefore, from (jHTj) . the Dirac eigenenergies corresponding to the bound-state 
solutions must be within the limits 

k^ 
E > -mc^ - 2 7^, for fci > (46) 

E < -mc^ + 2-^, for ki < (47) 
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For ki = 0, when there is a pure pseudoscalar potential hnear in x (the two-dimensional 
Dirac oscillator), Eq. (ji^ reduces to 

E^ = rr?c^ + he [(2n + 1) l^a] + fcs] • (48) 

Note that n > for A;2 > and ri > 1 for A;2 < 0, because for the latter case the lower 
component is proportional to a Hermite polynomial of degree n — 1. This result, together 
with Eq. (PS|). allows us to conclude that the Dirac eigenvalues are given hy E = ±Eq, where 



Eq = A/r?7?c^H-2(n+Ty^cpJ2| > rnc^, n = 0, 1, . . ., irrespective of the sign of k2. This means 
that the spectrum is independent of the sign of k2- The solutions E = mc^ and E = ~mc^ 
correspond to isolated solutions for /c2 < and ^2 > 0, respectively (see Eqs. (plj) and (jSH))- 
However, the spinors do depend on the sign of ^2, because different signs induce different 
node structures, as can be seen from P3|) . Since the charge-conjugation transformation 
changes the sign of ^2, it connects both kinds of solutions, but each can have both positive 
and negative energies. Thus a complete set of states for a certain k2 includes both particle 
(positive energy) and antiparticle (negative energy) states but they are not related to each 
other by a charge-conjugation transformation. The spectral gap between these solutions is 
2^m?c^ + 2hc\k2\ and does not disappear when m = 0. Therefore, there is no room for 
transitions from positive- to negative-energy solutions and Klein's paradox never comes into 
play. There are no scattering states when ki = 0, but for ki y^ the continuum states are 
those ones which have energies beyond the limits expressed by pF)|) - ()47j) . for which K < 0. 
As far as the chiral transformation is concerned, one can note that 7^ acts much in the same 
way as the charge- conjugation transformation, except that it does not change the sign of the 
energy and changes the mass sign. Therefore, the chiral transformation connects solutions 
with the same energy sign but that have opposite ^2 signs. 

When k2 = 0, i.e., no pseudoscalar potential, there is no isolated bounded solutions and 
the quantization condition expressed by Eq. ()44j) takes the form 

[E - mc^) a/(E + mc^) sgn(/i;i) = sgn(A;i) ( ^ + - ) hc^/2\kl\ (49) 

where sgn(/ci) stands for the sign function. This means that when ki > there are only 
bound states for fermions with E > mc^ and those bounded solutions are kept apart from the 
negative-energy continuum {E < —mc^) by an energy gap greater than 2mc^. On the other 
hand, for /ci < there are only bound states for antifermions with E < —mc^ and those 
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bounded solutions also do not mix with continuum states with E > —mc^. Therefore, the 
positive-energy bound states of fermions never sink into the Dirac sea of negative energies, 
and the negative-energy bound states of antifermions never reach the upper continuum, so 
again there is no danger of reaching the conditions for Klein's paradox. It is important to note 
that the solutions with fci < are not the antiparticles of the solutions with fci > because 
charge conjugation changes S into —A and vice-versa. This means that the antiparticles of 
the states with S = \kix^, A = 0, Ip = and /ci > are the negative-energy solutions of 
the relativistic system with A = — |A;ix^, S = 0, V^ = 0, with ki positive. 

To conclude, let us present, as an example of the rule stated at the beginning of the 
section, the eigenvalue equation for the relativistic harmonic oscillator when A = ^kix"^, S = 
0, V^ = k2x: 

E^ - m^c^ = (2fi + 1) he J ~^^ ki + kl-hek2 . (50) 

Here ft is the number of nodes (the degree of the Hermite polynomial) for the lower com- 
ponent, since the chiral transformation interchanges the upper and lower components, so 
that the second-order equation to solve in this case is the one for the lower component, 
Eq. (j29p . The upper component would be obtained from the corresponding first-order equa- 
tion Eq. (f^Hj) . As a final remark, one may note that for massless particles, the A = and 
S = oscillators have the same spectrum with the sign of k2 reversed. 

In 3+1 dimensions the corresponding eigenvalue equations to Eqs. (jl^ and (j^Uj) have a 
very similar structure. From [l^ we see that, for A = 0, one just adds the Dirac equation 
quantum number k to the coefficient of ^2 and modifies the coefficient of the square root to 
include the orbital angular momentum / of the upper component 



E"^ - Tr?e^ = {An + 2/ + 3) he J ^ ^^ ki + kj + {2k - l)hek2 , (51) 

while for S = one uses the orbital angular momentum / of the lower component 



E — mc^ 



^2 _ ^2^4 _ U^ ^ 21 + 3) he J ^^ki + kj- {2k - l)hek2 , (52) 



where k = —k (see |l5[)- If we compare these equations with ()44|) and ()50|). we see that 
the main features are identical, and thus all the main results obtained in the next section 
remain valid for the relativistic harmonic oscillator in three spatial dimensions. 
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IV. GRAPHICAL DETERMINATION OF THE ENERGY LEVELS 



We turn now to the solutions of Eq. ()44|) . In order to obtain a deeper insight on 
those solutions for arbitrary ki and /c2 one has to seek a convenient method for solving this 
transcendental equation. If we square Eq. (jiij) the resulting equation is in general a quartic 
algebraic equation in E, which in principle can be solved analytically. The price paid is that 
some spurious roots can appear in this process, although, of course, these can be eliminated 
by checking whether they satisfy the original equation. A more instructive procedure is 
follow a graphical method, by which one seeks the intersection points of the functions of 
energy in Eq. (j44|) : a parabola on the left-hand side 

fj,{E) =E^- m\^ , (53) 

and a square-root function on the right-hand side 



]s{E) = {2n + l)hc^{E + mc'^)ki/2 + kj + hck2 . (54) 

In this way, it is easy to see that there are at most two acceptable solutions of (|44|) . since 
the parabola and the square-root function can have at most two intersection points. The 
construction of the graph of fs{E) is split into three distinct classes depending on the ki 
values. Below we discuss these three classes of problems. It should not be forgotten, though, 
that the solutions with E = —mc^ are not included in this discussion, since they correspond 
to isolated solutions for ^2 > 0. 

This kind of analysis and the conclusions that will be drawn from it in the following 
subsections can be applied with appropriate modifications to the 3+1 case, and therefore to 
nucleon and antinucleon single-particle spectra. This is especially true as far the relations 
between the spectra of systems with A = and S = are concerned, since these are obtained 
from the general properties of charge-conjugation and 7^ transformations which hold in 3-1-1 
dimensions. 

A. The case ki = 

Although we have already solved this case, it is presented here for the sake of complete- 
ness. For ki = 0, fs{E) is just a nonnegative constant function as can be seen from Eq. 
having two symmetric intersections with the parabola. 
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Fig. 2] and Fig. |21 present the results for the first four quantum numbers. In these 
and the following plots we use the scaled quantities e = E/{mc^), ki = h'^ki/{m^c^) and 
K2 = hk2/{'m?c'). The first shows how the values of the energy are obtained from the 
intersection of the two curves and the second how the eigenvalues vary with k2. The energies 
for the zero mass case can also be computed from the intersections with the dashed parabola 
in Fig. ^ From these figures it is clear the symmetry of the energy levels. Note again that 
the isolated solution is not plotted. 
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FIG. 1: Scaled energies e = E/{m(P') for the first four levels when |k2| = 1 and ki = 0. The black 
dots represent the values of energies coming from the intersection of the parabola and the horizontal 
lines correspond to the scaled function fs{E)/{mc^)'^ for n = 0, 1, 2, 3. The dashed parabola is the 
function fp{E) for m = 0. In this case, the horizontal lines are just fs{E) for n = 0, 1, 2, 3, |k2| = 1 
in appropriate units and the horizontal scale is in the units of the energy E. 



B. The case /ci > 



When ki > 0, fs{E) is a monotonically increasing function of E. For each n, it can 
intersect the parabola at one or two points, depending on the particular values of fci and k2- 
To determine the conditions for this to occur, we have first to determine the minimum of 
fsiE), which occurs when the expression under the square root in ()54|) is zero: 

h 



E^ 



-mc 



2^. 



(55) 
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FIG. 2: Scaled energies e = E/imc?) for the first four levels when ki = as a function of K2. For 
m = the plots are qualitatively the same, except that the lines for the different levels begin at 
the origin and the vertical scale would be in units of the energy E. 



This value sets a lower bound for the energies of all the solutions of Eq. 
the function fs{E) takes the value 



fs ■ = hcko 



For E = Emin, 



(56) 



and /„ takes the value 



J Prain 



, rCo / nJ 



+ mc 



(57) 



ki \ki 

When ki > 0, the minimum of fs{E) occurs for negative values of the energy, that is, on 
the left (and decreasing) side of the parabola. If fp^.^^ < fs^^^ there is only one solution 
of Eq. P^ and it has positive energy. From ()56|) . this can happen only for k2 > 0, since 
fpmin > 0- ■'■f /smin — /pmin there are always two solutions, one negative and one positive. 



because /s(0) = {2n + l)hc\Jmc'^ki/2 + /i;| + hck2 is always positive, no matter how large 
is I/C2I, for a positive ki. Furthermore, the positive and negative solutions are such that 
\E\ > mc^, since, when \E\ < mc^, fs{E) > fp{E) and thus there are no solutions. Note, 
however, that for n = and k2 negative there is a solution with E = —mc^, but this is an 
isolated solution out of the Sturm-Liouville problem we are considering. 

The considerations of the previous paragraph can be further elaborated by analyzing the 
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behavior of the function F{k2) 

F{h) = U_ - /._ = ^ (kl + mc^hk, - ^) . (58) 

According to what was stated in the previous paragraph, if F{k2) > one has two solutions 
for each n, and if F{k2) < there is just one solution for each n. To determine the signs of 
F{k2) for each value of A;2; we will need to know the values of k2 which are roots of F{k2) = 0, 
or of the equation 

kl + mc^kik2 -^ = , (59) 

excluding the solution A;2 = 0. Applying Descartes' rule of signs to this last equation, we see 
that it can only have a positive real solution. The rule states that the number of positive 
real roots rtj. of an algebraic equation with real coefficients 

auz^ H h aiz + ao = (60) 

is never greater than the number of changes of signs ric in the sequence Ofc, . . . , Oi, Oq (not 
counting the null coefficients) and, if less, then differs from ric by an even number. This 
means that rtr can take the values n,, = "^.c, nc — 2, nc — 4, . . . , Uc — 2n > 0, n = 0, 1, . . . . 
In this case, for positive fci, there is only one such sign change and therefore there is only 
a positive and unique solution, which we shall call k"^. The function F{k2) changes sign at 
/c2 = and ^2 = k2, being negative in the interval < k2 < k2 and positive outside this 
interval. Therefore, according to what we have said before, there is only one positive-energy 
solution when < A;2 < ^2 ^^^ ^wo solutions (one positive and one negative) when k2 < 
or k2 > k2- li k2 = 0, the two functions intersect at -Emin = —mc^, but again, since this 
solution was excluded from the Sturm-Liouville problem leading to the eigenvalue equation 
()44|1 . we consider only the harmonic oscillator solutions with energy greater than mc^. 

The graphical representations of fp and /^ are shown in Fig. 01 for m ^ 0. We plot the 
square- root function for three values of K2 (the scaled ^2), for a positive ki. We plot fs{E) 
with n = for K2 = 0, 0.6 and with n = and n = 1 for ^2 = —0.6. The value of Ki 
was chosen such that -^(^2) has a positive value for K2 = 0.6. There is always an infinite 
and unbounded set of positive-energy states with E > mc^ for any value of /c2, and there 
are always negative states when k2 is negative. Note that increasing the value of n does 
not change qualitatively these features, as can be seen from the plots of fs{E) when ^2 is 
negative. 

21 



\ 




8 






A- 


— K^ = 0.6 

-K^= 

.- K =-0.6 


6_ 


_/■ / 
^ ' / 

y' / 

.^ n=l / 

y j 




\ 






y / 








\ 




y / 




\ 


*y / 






\ / 


/ 


. — ■ 




1 .ir^^'^o / 1^. 


" n=0 




-2 ^X_ 


y\ ' ▼" w ^ 


-3 


__^^/ 1 2 3 






e 



FIG. 3: Graphical solution of equation H44() for ki = 4 and three values of K2- The circle, square 
and diamond symbols mark the energy solutions for n = and K2 = 0.6, 0, —0.6, respectively. In 
the case of K2 = —0.6, the square root function for n = 1 is also represented and the respective 
energy solutions are marked by the triangle symbols. 

In Fig.|3]it is clearly seen that for values of k2 between zero and fcg there are only positive- 
energy solutions. Moreover, one can see from Eq. (jSUj) that the value of fcg increases with fci, 
meaning that when ki ^ \k2\ one reaches the limit of the harmonic oscillator (fci > 0, A;2 = 0) 
where there are only positive bound-state solutions with E > mc^. Thus, the existence of 
the pseudoscalar potential gives rise to negative-energy solutions coming from the Dirac sea. 
For A;2 < these kind of solutions always exist. We excluded the isolated level with n = for 
negative ^2. Note also that for decreasing negative /c2 the positive-energy solution for n = 
goes to E = mc^, which was the case we discarded in the discussion of the Dirac oscillator 
in the previous section. When ^2 = ^2 there is an infinite set of degenerate solutions with 
E = -Emin, because in this case the intersection point of fs{E) with the negative-energy 
branch of fp{E) is the initial (minimum) point of the square-root function, which of course 
does not change with n (see Fig. 121). However, when E = -Emin, the value of K defined in (IH7|) 
is zero, and thus, from ^^, A = and therefore these solutions do not correspond to bound 
states (see Eqs. fj39|l and (jl2|) ). For ^2 ^ ^2 the degeneracy is removed and we have a very 
high density of very delocalized states (since A > 0) which then decreases with increasing /c2. 
In Fig. m we also plotted the parabola —mc^ — 2/i:|//ci (in scaled units —1 — 2/i:2/fi;i) which 
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corresponds to the asymptotic values of the negative-energy states when n -^ oo and thus 
represents the set of lower energy bounds of the negative-energy solutions for each pair of 
ki and ^2 values. This parabola can be also seen as delimiting the region in the energy-fc2 
plane in which K <0, i.e., where there are no bound-state solutions. This is related to the 
discussion in the previous section leading to Eqs. (|iH|) and (jTTj) . 

As referred before, for a given energy E < ~mc^ and positive ki there are not bound 
states for values of ^2 in the interval —\J—{E + mc^)ki/2 < k2 < \/—{E + mc^)ki/2. In 
fact, from what was said above in this section, one can be more precise and state that there 
are no negative bound states for < /c2 < A;2- Since /cg is zero when fci = and increases with 
increasing /ci, this means that when fci = (the Dirac oscillator) there are no continuum 
states and when ki -^ oo (or when A;2 "~* 0), all the negative-energy states are continuum 
states. 




FIG. 4: Scaled energies for the first six levels when ki = 1 as a function of K2. The n = level 
for negative /c2 and negative energy was not plotted because it corresponds to the isolated solution 
with e = —1 {E = —m(?). The value of K^ = hk2/{m?(?) is marked with a cross in the energy 
axis. The hollow circles mark the points {k2 = 0, e = — 1) and {k2 = k2, e = Cmin = E^i^/mcP' for 
which there are no bounded solutions. The dashed parabola represents the function —1 — 2^2/1^1, 
whose meaning is given in the text. 
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C. The case ki < 

For ki < the function fs{E) decreases monotonically with the energy, reaching again 
its minimum when the expression under the square root in (j54j) is zero, setting now a 
maximum value for the energy of the solutions of Eq. (j44|) . The expressions for -Emax, 
/smin = /s(-Emax) and fp^.^ = fp{E^s.yi) are the same as corresponding ones in the previous 
subsection, respectively Eqs. (j33j) . (j3Bj) and (jSTj), except that, of course, now ki is negative. 

Noting that fs ai E = —mc^ has a nonnegative value given by fs{—mc^) = {2n + 
l)hc\k2\ + k2, we see that there is always an infinite and unbounded set of negative-energy 
solutions (one for each n) with E < —mc^. For n = and k2 < 0, there is a solution with 
E = —mc^ which was already discussed in the previous section. However, when ki < 0, 
one can solve Eq. ()30|) when E = —mc^ and obtain a bound-state solution with both upper 
and lower components different from zero, if and only if k,2 = k,i/8. One can show that 
this corresponds to a ra = solution given by the functions ()39|) and (|1^ taking the limit 
E -^ —mc^. 

In addition to that set of solutions, there can be another set of solutions only if the 
function E{k2) (see Eq. I3S|) is zero or positive, as explained before. To determine the values 
of /c2 for which this happens, one has to compute the roots fcg of E{k2) = 0. When ki < 0, 
Descartes' rule of signs applied to E{k2) = states that one can have a positive root and 
two or zero negative roots (this is done by changing z -^ —z in ()60|) ). in addition to, of 
course, the root k2 = 0. The two negative roots arise for a value of ki = ki such that 
-^'(^2) ~ 0- ■'■^ scaled units, this value is kI = —64/27, corresponding to the double negative 
root ^2 = —8/9. For Hi < k,\ there is only a positive root, while for k\ < ki < {] there are 
two negative roots and a positive root. Because -^(^2) changes sign at k2 = ^2 (except when 
^2 is a double root) and E{k2) is always positive when |/c2| —* 00, this means that when 
Ki < kI there is only one range of values of k2 where -^(^2) < whereas when ni > kI there 
are two ranges of ^2 values where this happens. Within these ranges, only negative-energy 
solutions exist. When ki = nl, the value ^2 = ^^2 = —8/9 corresponds to an infinite set of 
degenerate negative-energy solutions with scaled energy e = — 1 — 2{k,2)'^/k,'1 = —5/3, which, 
as mentioned before, are not bound-state solutions. These features of the spectrum can be 
seen in Figs. E] and El 

In Fig. El the solutions with —mc^ < E < with n = 1, . . . , 5 are basically degenerate 
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FIG. 5: Scaled energies for the first six levels when ki = — 1 > kJ as a function of K2- As in Fig. 0] 
the three non-zero values of the roots of F(k2), K2, are marked with a cross in the energy axis. 
As in that figure, the hollow circles mark the point (k2 = 0, e = — 1) and the three points where 
(A;2 = A;|, e = Cmin)- The dashed parabola represents the function —1 — 2^2 //^i (see text). 




FIG. 6: Scaled energies for the first six levels when ki = —3 < k^ as a function of k.2- The 
conventions and meaning of the symbols are the same as in Figs. |3and[5J The dashed parabola 
represents the same function as in Fig. |SJ 

(they follow closely the parabola —1 — 2/T;2/fi;i) and the n = solution is very close to them. 
Since, as remarked before, on the parabola the value of K is zero, these are very weakly 
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bounded solutions. Note that the n = positive-energy solution with E < mc^ goes to mc^ 
when k2 -^ — oo, thus becoming the isolated solution of the Dirac oscillator {ki = 0). This 
can also be seen in Fig. IHl where one may also remark that there is a value of K2 around —0.9 
for which the solutions for distinct values of n and —mc^ < E < are almost degenerate. 
This point corresponds to the local minimum of E{k2), which, for Hi = —3, is very close to 
the double root Kg = —8/9. Similarly as in the case of ki > 0, the solutions with A;2 = ^2 7^ 
have infinite degeneracy, but they do not correspond to bounded solutions. Again, for k2 
near fcg there is a high density of very delocahzed states. The parabola —1 — 2K2I ki sets 
now an upper bound {n -^ 00) for the positive-energy states. 

The limits /cg — * and \k2\ -^ 00 correspond, respectively, to the harmonic oscillator 
solutions for negative-energy (antiparticle) states and to the Dirac oscillator. In the first 
case there are no bound states for particles (positive-energy states). As \k2\ increases, 
positive bound states start to appear. 

D. Variation with ki for /c2 fixed 

To complete the study of solutions of equation (jl^ we present now a plot of its solutions 
when ^2 is kept fixed and let ki vary. Note that the existence of one or two solutions is 
governed by the sign of the function 

G{k,) = ^ {kl + m^k,k2 - ^) , (61) 

identical with F{k2), except that now k2 acts as a parameter and ki is the variable. The 
critical value is now Kg = — 1. In Fig. Owe present a plot of the solutions for /c2 > 0, when 
there are two possible values for kl, the roots of G{ki) = 0. Apart the factor 4k2/kl, this 
is just a quadratic equation in ki, so that it is easily shown that when /t:2 < — 1 it has no 
real solutions. This means that, since in this case G{ki) > 0, there are always two sets of 
solutions. Note that now the points in the energy-fci plane for which K = are on the 
hyperbole —1 — 2fc|/fci, such that only in the region between the hyperbole branches, where 
K is positive, bound states may exist. 

As before, we can recognize the Dirac oscillator and the harmonic oscillator limits, as one 
sets kl = and |A;i| — *> 00, respectively. In this last case, one has the harmonic oscillator 
solutions for particles or antiparticles depending whether ki goes to +00 or —00. One can 
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FIG. 7: Scaled energies for the first six levels when K2 = 1 > /i^i ^s a function of ki. As before, the 
hollow circles mark the energy solutions for ki = /cjj' . The dashed hyperbole represents the function 
— 1 — 2^2 /ki (see text). 

again see that the energy levels of the Dirac oscillator are symmetric about E = and 
that there are no scattering states. Here this follows simply from the fact that A;i = 
is an asymptote of the hyperbole branches which separate the continuum region from the 
bound-state region. 

V. SOLUTIONS FROM SYMMETRY ARGUMENTS 

As referred before, the case S = can be obtained from the A = case by applying 
the charge-conjugation transformation. We recall that this transformation performs the 
changes E — > — E', A -^ — S and Vp -^ —Vp, or, for harmonic oscillator potentials, the 
changes ki -^ —ki and k2 -^ — /c2- This can be seen if we solve the eigenvalue equation for 
S = 0, Eq. (|^. for ki = 1 and plot the solutions for several values of K2. The result is 
shown in Fig. |H1 

If we compare this figure with Fig. El (/ti = —1) we see that the plots are identical if we 
reverse both the vertical and horizontal axes, i. e, if we reverse the sign of the energy and 
of ^2, respectively. Of course, the identification of the particle and antiparticle levels is also 
reversed, as it should be, since we are applying the charge-conjugation transformation. 
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FIG. 8: Scaled energies for the first six levels when ki = 1 as a function of H2 and S = 0. The 
conventions are the same as in Fig.|Sl The dashed parabola represents now the function 1 — 2^2/^1. 
n stands for the quantum number of the solutions of Eq. I|50)l. 

The existence of both positive- and negative-energy solutions in a system with S = can 
be relevant for nuclei. In the 3+1 Dirac equation of relativistic nuclear mean-field theories 
there is a connection between the (isoscalax) vector and tensor potential. The latter is 
proportional to the derivative of the former [23| • For the harmonic oscillator potentials used 
in this paper, this would give, using the notation of Ref. [2^1 , the relation 



fi2 = ^1^1 ■ 



(62) 



The constant f^ can take values up to around 1.3 in nuclei using relativistic mean-field 
potentials of Woods-Saxon type 



|. For harmonic oscillator mean- field potentials, there 
should be also an upper value for /t,, and therefore the relation ()fi2|l sets a maximum for K2 
in nuclei. As was seen above, this is relevant to know whether there can be simultaneously 
nucleon and antinucleon bound states in nuclei. This would depend very much on the 
strength of the A potential, which is given by ki. As is apparent from Fig. |H1 for small 
values of the scaled fci, ki, the maximum value of K2 allowed by Eq. (j62j) could be such that 
there would exist only positive-energy states. On the other hand, for stronger A potentials, 
the relation (jH^ could allow both positive- and negative-energy bound states, as can be seen 
by applying the axis reversal described above to the A = spectrum shown in Fig. IHl 
In Sec. mil we have seen that the 7^ chiral transformation can also be used to get formally 
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E = from A = solutions, involving, however, the unphysical transformation m -^ —m. 
This objection is overcome when we deal with ultrarelativistic particles, for which we may 
consider m ~ 0. In this case, the chiral transformation relates massless physical states with 
A = 0, S = y and a pseudoscalar potential Vp to massless states with S = 0, A = l^ 
and —Vp, where V can be an arbitrary potential. For harmonic oscillator potentials, one 
has E = kix'^/2 ^ A = kix'^/2 and k2 — > —k2. Since, as explained in Sec. IIIH the chiral 
transformation does not change the equations of motion, the energy of their solutions is not 
changed. This can be seen in Figs. Inland ^J where are plotted, respectively, the solutions of 
the A = and S = eigenvalue equations, Eqs. ()44|1 and ()5()j) . when m = 0, ioi k, = 1 and 
several k2 values. The basic features of these plots were already discussed in the previous 
section. As before, the E = isolated solution was excluded. We see that the spectra are 




FIG. 9: Energies for the first six levels when ki = 1 as a function of K2 for A = and m = 0. The 
conventions are the same as in Fig. [l] The dashed parabola represents the function — 2k|/k;i. 

the same apart from the k2 -^ — ^2 change, which shows up as a reversal of the horizontal 
axis. For ^2 = 0, the pure harmonic oscillator solution, both spectra are identical, with no 
negative-energy solutions and, of course, no nonrelativistic limit. Note, however, that the 
A = and S = energy eigenstates are not the same, since 7^ interchanges the upper and 
lower components. This amounts to say that these states are not eigenstates of 7^, which of 
course they could not be, since, because oi Vg, it does not commute with the Hamiltonian 
((Tj) even with m = and Vp = 0. If Vs would be zero, then, because one has either A = 
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FIG. 10: Energies for the first six levels when /ti = 1 as a function of K2 for S = and m = 0. 
The conventions are the same as in Fig. 01 The dashed parabola represents the function — 2K2/K1. 

or S = 0, Vt would have also to be zero, and we would have a free ultrarelativistic particle, 
which has a good chirality, i.e., is an eigenstate of 7^. 

VI. CONCLUSIONS 



In this work we studied the solutions of 1+1 Dirac equation with a potential with the 
most general Lorentz structure which leads to harmonic oscillator potentials in the second- 
order differential equation for the spinor components. This is achieved when the vector Vt 
and scalar Vs potentials are quadratic and satisfy either A = Vt — Vs = or T, = Vt + Vs = 0, 
and the pseudoscalar potential is linear. We showed that there are no bound-state solutions 
in the nonrelativistic limit for E = and Vp very small, meaning in this case we have an 
intrinsic relativistic problem. 

We analyzed in great detail, by use of a graphical procedure, the solutions of the tran- 
scendental eigenvalue equation in the case of A = 0, for all signs of the potentials and for 
both particle and antiparticle states. As referred in the Introduction, there is an increas- 
ing interest in studying the behavior of antinucleon in nuclei, and in particular how well 
pseudospin and spin symmetries hold for nucleons and antinucleons respectively. 

We discussed in detail the isolated solutions, i.e., solutions out of the Sturm-Liouville 
problem, as well as the critical points in the parameter space of ki and k2, the coefficients 
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of the quadratic and linear potentials. Near these points there is a high density of bounded 
states which are very extended in space. These critical points are also relevant to determine 
which values of ki and k2 allow to have a spectrum of both fermion and antifermion bounded 
solutions simultaneously or just one of these kind of solutions. We established a connection to 
3+1 dimensional relativistic oscillator harmonic potentials applied to nuclei, where there is a 
relation between those coefficients and concluded that, for a sufficiently strong A potential, 
both kinds of solutions {i.e, for both nucleons and antinucleons) can coexist. 

We showed how to obtain the solutions for S = from the A = case, using the charge- 
conjugation and chiral transformations. In particular, we used the chiral transformation to 
show that, for massless fermions and zero pseudoscalar potential, the A = and S = 
solutions have the same spectrum. These conclusions are quite general and do not depend 
on the number of dimensions so that they remain true in 3+1 dimensions. Also, most of the 
features of the 1+1 spectra presented in this work are also present in the 3+1 relativistic 
harmonic oscillator. Indeed, since the 1+1 eigenvalue equations for A = and S = are 
very similar to the corresponding 3+1 equations, we can apply the same graphical methods 
to the 3+1 case, and therefore most of the analysis of the 1+1 case remains true with 
appropriate changes. For instance, the shape of the parabola in the energy-A;2 plane related 
to -Emin would depend now on the value of the spin-orbit quantum number k. The function 
F{k2) would also depend on k, and therefore the roots fcg and critical value fcg would also 
depend on k. For A;2 = there would be the known harmonic oscillator 2n + 1 degeneracy, 
produced by the U(3) symmetry studied very recently by Ginocchio [22| in connection to 
spin and isospin symmetry. 

Finally, the chiral transformation enables us to switch between the so-called pseudospin 
(S = and Vp = ^i) and spin symmetries (A = and ^ = 0) ja |l^ for any potential. 
In the context of a study of the symmetries of a 3-dimensional relativistic harmonic oscil- 
lator, Ginocchio already showed that one can transform the generators of the SU(2) spin 
symmetry into the pseudospin generators by a 7^ transformation J22|- We showed that this 
transformation allows us to relate the A = 0, S = V and S = 0, A = V problems for any 
potential V , and, for massless fermions and zero pseudoscalar potential, also to conclude 
that the two energy spectra are the same. 



31 



Ackno^vledgments 

We acknowledge financial support from CAPES (PDEE), CNPq, FAPESP, and FCT 
(POCTl) scientific program. 



[1 
[2 
[3 

K 

[5 

[6] 
[7] 
[8 
[9 

[lo; 
[11 

[12 
[13 
[14 
[15 

[16 

[17 
[18 
[19 

[2o; 



T.-S. Chen, H.-F. Lii, J. Meng, S.-Q. Zhang, and S.-G. Zhou, Chin. Phys. Lett. 20, 358 (2003). 

V. I. Kukuhn, G. Loyola, and M. Moshinsky, Phys. Lett. A 158, 19 (1991). 

J. N. Ginocchio, Phys. Rev. C 69, 034318 (2004). 

J. N. Ginocchio, Phys. Rev. Lett. 78, 436 (1997). 

J. Meng, K. Sugawara-Tanabe, S. Yaniaji, P. Ring, and A. Arinia, Phys. Rev. C58, R628 

(1998). 

J. N. Ginocchio, Phys. Rep. 315, 231 (1999). 

S. Marcos, L. N. Savushkin, M. Lopez-Quelle, and P. Ring, Phys. Rev. C 62, 054309 (2000). 

J. N. Ginocchio, Phys. Rep. 414 165 (2005). 

J. Meng, K. Sugawara-Tanabe, S. Yaniaji, and A. Arima, Phys. Rev. C 59, 154 (1999). 

P. Alberto, M. Fiolhais, M. Malheiro, A. Delfino, and M. Chiapparini, Phys. Rev. Lett. 86, 

5015 (2001). 

P. Alberto, M. Fiolhais, M. Malheiro, A. Delfino, and M. Chiapparini, Phys. Rev. C 65, 

034307 (2002). 

R. Lisboa, M. Malheiro, and P. Alberto, Phys. Rev. C 67, 054305 (2003). 

R. Lisboa, M. Malheiro, and P. Alberto, Braz. J. Phys. 34, 293 (2004). 

J. S. Bell and H. Ruegg, Nucl. Phys. B98, 151 (1975). 

R. Lisboa, M. Malheiro, A. S. de Castro, P. Alberto, and M. Fiolhais, Phys. Rev. C 69, 024319 

(2004). 

D. ltd, K. Mori, and E. Carriere, Nuovo Cimento A 51, 1119 (1967). 

M. Moshinsky and A. Szczepaniak, J. Phys. A 22, L817 (1989). 

A. S. de Castro, Ann. Phys. (N.Y.) 316, 414 (2005). 

S.-G. Zhou, J. Meng, and P. Ring, Phys. Rev. Lett. 91, 262501 (2003). 

T. D. Cohen, R. J. Furnstahl, and D. K. Griegel, Phys. Rev. Lett. 67, 961 (1991). 



32 



[21] P. Strange, Relativistic Quantum Mechanics with Applications in Condensed Matter and 
Atomic Physics (Cambridge University Press, Cambridge, 1998). 



[22; 
[23 

[24 
[25; 
[26 
[27; 
[28 
[29 
[30 
[31 
[32; 
[33 
[34 
[35; 
[36 

[37; 

[38 
[39 



J.-Y. Guo, X.-Z. Fang and F.-X. Xu, Nucl. Phys. A 757, 411 (2005). 

R. Lisboa, M. Malheiro, A. S. de Castro, P. Alberto, and F. Fiolhais, Int. J. Mod. Phys. D 

13, 1447 (2004). 

Y. Nogami and F. M. Toyama, Can. J. Phys. 74, 114 (1996). 

F. Toyama, Y. Nogami, and F. A. B. Coutinho, J. Phys. A 30, 2585 (1997). 

F. M. Toyama and Y. Nogami, Phys. Rev. A 59, 1056 (1999). 

R. Szmytkowski and M. Gruchowski, J. Phys. A 34, 4991 (2001). 

M. H. Pacheco, R. Landim, and C. A. S. Almeida, Phys. Lett. A 311, 93 (2003). 

A. S. de Castro, Phys. Lett. A 318, 40 (2003). 

A. S. de Castro, Ann. Phys. (N.Y.) 311, 170 (2004). 

B. Thaller, The Dirac Equation (Springer-Verlag, Berlin, 1992). 

C. Itzykson and J.-B. Zuber, Quantum. Field Theory, (McGraw-Hill, New York, 1980). 
A.S. de Castro, Ann. Phys. (N.Y.) 320, 56 (2005). 



A.S. de Castro, |hep-th/05 10003] Int. J. Mod. Phys. A, in press. 



L. I. Schiff, Quantum Mechanics (McGraw-Hill, New York, 1968) 3rd ed. 

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Dover, Toronto, 

1965). 

I. N. Bronshtein and K.A. Semendyayev, Handbook of Mathematics (Springer-Verlag, Berlin, 

1998). 

P. Alberto, R. Lisboa, M. Malheiro, and A. S. de Castro, Phys. Rev. C 71, 034313 (2005). 

J. N. Ginocchio, Phys. Rev. Lett. 95, 252501 (2005). 



33 



